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We discuss general formation of complementary behaviors, functions and forms in 
biological species competing for resources. We call orthogonalization the related 
processes on macro and micro-level of a self-organized formation of correlations in 
the species properties. Orthogonalization processes could be, for example, easily 
observed in sympatric speciation, as we show in numerical studies carried with a 
new population equation. As a practical result, we find that the number of species 
is proportional to the effective richness of resources and depends on their history. 



Introduction 

Individuals in any biological species differ a little in their behaviors, forms and 
other parameters, which can help to relax their otherwise large competition. If the 
differences become too large, individuals can decrease their interest for interbreeding, 
and species can consequently split up. During speciation and species coexistence 
correlations are formed between parameters of the species, which adjust distributions 
for different species in such a way that resources are effectively shared. In this process 
complex patterns of complementary sets of parameters are formed, where species can 
be seen to approach discrete orthogonal states. These "orthogonalization processes" 
(OP) can be modeled in multidimensional mathematical spaces, with coordinates 
given by the species parameters. 

It appears natural to define species by a mutual reproducibility of their indi- 
viduals (Mayr (1942, 1963)). Since OP are strong in clusters of evolutionary and 
functionally close species, they could give a clear substance to the species defini- 
tions. Globally, we can define species orthogonalization as a process in which species 
approach a state, where they share resources and habitats in the most effective way, 
fully stabilized in time. Since species not only passively approach this ideal state, 
but they prefer to overturn the "status quo" , by permanently improving the quality 
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of their gene equipment, full orthogonalization can not be realized. Instead, evolu- 
tion gives a relatively stable coexistence (quasi-orthogonality) of biological species, 
time from time disturbed by speciation events. Close to these events OP are the 
most intense, but they are strong until the new species "settle down". Species in- 
teractions, which enable this process, last even after the species get in the resulting 
orthogonal states. 

The need for orthogonalization might have caused the disappearance of many 
"intermediate" species, since these would hamper orthogonalization in the short 
and long term run, the last because both the resources volume and speciation rate 
are relatively stable. A more specific trace of OP are character shifts observed in 
species sharing habitats or in sympatric speciation (Brown & Wilson (1956), Smith 
(1966), May & Mac Arthur (1972), Slatkin (1980), Schluter & Mac Phail (1988), 
Doebeli (1996a), Drossel & Mac Kane (1999)), which reveal efficient division of 
resources from the point of structure, day time or form of consumption. Competition 
for (orthogonal) discrete niches (Hutchinson (1968)) is confirmed, for example, by 
simple speciation/extinction models (Valentine & Walker (1988)). Therefore, OP 
might even control the numbers of species living in given areas (Rosenzweig (1995), 
Plotkin et al. (2000)). 

To understand OP, we discuss first some evolution aspects of species on the 
micro and macro-level, and flow of information between the levels. Then, we focus 
on modeling of a sympatric speciation of species (Doebeli (1996b), Bagnoli & Bezzi 
(1997), Kondrashov & Kondrashov (1999), Dieckman & Doebeli (1999), Drossel & 
Mac Kane (2000)), and find a new speciation equation. We present it in numerical 
examples, where OP can be observed. 

The concept of species orthogonalization 

We can look on species orthogonalization from different levels of the OP activity. In 
a self-organized molecular system (Eigen (1971)), different types of molecules play 
different roles, expressed by their microscopic structures. Evolution processes lead- 
ing to separation of these roles can be called orthogonalization on the micro-level 
(OMIL). Since correlations are also built between clusters of cooperating but dif- 
ferent types of molecules, natural orthogonalizing units are molecular quasispecies 
(Eigen & Schuster (1977), Eigen, M., McCaskill, J, k Schuster, P., (1989)). There- 
fore, eigensolutions of the Eigen's equations implicitly incorporate quasispecies or- 
thogonalization on the micro-level. 

The molecular system can be a part of a macro-system, competing or coop- 
erating with other macro-systems of the same or different types. Different types of 
macro-systems (individuals in different species or different organs inside individuals) 
analogously develop complementary functions and signaling, activities in different 
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times or space regions. Processes leading to such a differentiation of biological 
species can be called orthogonalization on the macro-level (OMAL). We could again 
find that OMAL is rather realized between self-organized clusters of species, anal- 
ogously to the quasispecies on the micro-level, which can help to improve species 
definitions. Forces leading to OMAL have predominantly macroscopic origin, but 
their influence must be "inheritably" fixed on a micro-level, which can be realized 
there as a part of "externally driven" OMIL. 

REALIZATION OF ORTHOGONALIZATION PROCESSES 

Realization of OP on all structural levels of biological species is rather fascinating. 
Starting from the micro-level, proteins on the tertiary level, for example, resemble 
macro-tools, where prototypical sections can be identified with letters in a hypo- 
thetical alphabet (Lesk (1951), Creighton (1992), Holm & Sander (1998)), showing 
the activity of OMIL, i.e. a functional orthogonalization of certain protein sections. 
Recently, for example, orthogonalization as a measure of similarity between differ- 
ent molecules was used in a modeling of prebiotic species formation (Segre et al. 
(2000)). Proteins also reflect the environment temperature, pressure or acidity by 
the compactness of their folding (Lumsden et al. (1997)), which shows that some de- 
grees of freedom can be tuned in proteins from the macro-level, without qualitatively 
affecting their functions. 

Comparative studies show (Creighton (1992)) that proteins, common among 
different species, have nearly zero mutations at some regions, reflecting functional 
importance of these sections. Other sections are nearly neutral with respect to 
protein's functions, since compositions of amino-acids present here are reasonably 
varied between evolutionary distant species (Wang (1996)). Such a structural "self- 
averaging" often applies to elongated sections stretching out of the proteins. The 
fact that fluctuations in these sections between individuals in a single species are 
still statistically small could reveal their use for ongoing OMAL processes. These 
processes can be engaged in a structural tuning of the neutral sections, with the 
goal to influence speed (catalytic strength), timing and other aspects of protein's 
functions, and thus vary the species morphology and general activity. Some genes, 
present in the gene-pool in many similar copies (Creighton (1992)), are also tuned 
by OMAL to perform the same functions in different evolution stages of the organ- 
ism, which carries them. Therefore, proteins and other microscopic units probably 
first internally develop by OMIL, and later are being externally tuned by OMAL, 
according to the needs on the macro-level. 

OMAL is realized through an information (vertical) flow between the species 
macro and micro-levels, where each level expresses its needs and possibilities in 
a specific way. Stimulations between the levels are rich: the macro-world provides 
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species individuals with various resources, acts on them through a set of macroscopic 
variables $ (like radiation, temperature, humidity, ...), and checks their ability to 
fit a given time-dependent environment. The individuals react according to their 
inherited and gained informations (Irwin & Price (1999)), coded in their gene-pools 
and memories, evolving as a response to the living conditions. 



Mathematically, an individual can be represented by a point in a vector space, where 
coordinates 0, are related to codons, genes, or other microscopic units (Baake et al. 
(1997)). Giving the fluctuations in their composition, a species with many indi- 
viduals occupies a finite volume A0 of this space. An extreme example are viral 
quasispecies (Eigen (1993)), where A<f) is very large and rapidly grows in time by 
mutations. The species individuals can be also described in a vector space with coor- 
dinates given by their macroscopic characteristics x (Priigel-Bennett (1997)). Some 
sets of microscopic parameters <fi (genotypes) can lead to morphological solutions 
with close macroscopic parameters x (phenotypes). Ambiguity of this reflection can 
be traced by a genotype-phenotype mapping (Fontana & Schuster (1998)). 

The two types of vector spaces are formed by a direct product of their subspaces 
7i = 7ii x 7i 2 x ... x 7i n for the involved microscopic or macroscopic parameters. 
In dependence on the studied level, an individual can be thus identified with the 
vector x = (xi, x 2 , x n ) or = (0i, 2 , 4>n)- A species with many individuals 
can be described by a (character) distribution /(x) or /(</>), which reflect various 
correlations in the populations. The distributions can be smoothed and normalized 
/o(x) = f(x)/N, where N is the number of involved individuals. 



To quantitatively appreciate OP, we can formally introduce a measure of species 
orthogonality. We say that two species are orthogonal with a weight /i, if a scalar 
product (Prugovecki (1981)) of their distributions (/o|<7o) fulfills 



Any species is orthogonal to itself with a weight \l = 0. Similarly, \x = 1 for any 
two species, if their distributions for at least one component of the vector x totally 
avoid each other. Definition (H) just touches the aspect of a "visual dissimilarity" 
of the character distributions in species. The value /x ~ 1, for example, does not 
assure that the two involved species are related in any way, and if they are, it does 
not prove that OP are completed. The last should be independently checked from 
the time-independence of the distributions /, g, as outlined in the Introduction. 



Description of orthogonalization processes 



EVALUATION OF ORTHOGONALIZATION PROCESSES 



(fo\go) = J 
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To get more refined tools for capturing OP, we should rather focus on some 
"dynamical aspects" of OP, like the fact that species living in close contacts often 
largely interact, even if they are practically orthogonal (f\g) ~ 0. Interactions, 
which enable orthogonalization, are realized through competition, described here by 
broadened distributions fc, gc ( see Eqn. g). For interacting species, these functions 
thus have a nonzero overlap (fc\gc) 7^ 0, even if (f\g) ~ 0. When the interactions 
proceed vertically on the food-web (Amaral & Meyer (1999)), OP can be traced 
analogously. 

It is also important to know in which way OP form correlations between the 
parameters xi, x 2 , ■ x n of interacting species. This can be tested from the species 
distributions F(xk,xi, ...), G(xk,xi, ...), resulting by projection of /(x), g(x) on the 
selected Xk,xi, ... parameters 

F, G(xk, xi, ...) = J dx\dx2...dxk-idxk + \...dxi-idxi + \...dx n f,g(x) ■ (2) 

Similarly as fc, gc, the projected distributions F, G could also overlap 

(F\G) = J dx k dx t ... \jF(x k ,x h ...)G(x k ,x h ...) ^ , (3) 

even if the species are quasi-orthogonal (f\g) ~ 0. For interacting species, this 
can be largely due to correlations formed by OP between the individual parameters 
x i,2,...,m m contrast to the overlap of fc, gc due to broadening (interaction). 

These correlations can be traced, if we project /, g on a functionally relevant 
cluster of parameters. An example for two interacting bird species might be the 
triplet: x^-beak length, x/-leg length and x^-wing size. This triplet could give a 
reasonably small overlap (F\G) ~ 0, while projections on subsets of (xb,xi,x w ) 
might already result in large overlaps (F\G) ^ 0. This can be especially relevant for 
the projected and competition broadened distributions F c , Gc- We can thus define 
an effective orthogonalization variable x, which would give the least overlap when 
we project on it. This variable then speaks about the character of correlations built 
by OP. The simplest possibility how to define it (see Fig. |l|) is to project along the 
vector x = (Cf — Cg) x, oriented in the direction between the centers of gravity 
Cf = {x^,xf,x^) and Cg = ( x b ■> x ? ■> x w) °f the projected species distributions 
F(xb,xi,x w ), G(xb,xi,x w ), and integrate these distributions out in the remaining 
two orthogonal directions. We will apply the effective one dimensional variable x in 
our numerical studies. 

In Fig. [I] we schematically present the distributions f(xi,x 2 ), g(xi,x 2 ) for 
two species. They are practically orthogonal (f(xi,x 2 )\g(xi, x 2 )) ~ 0, while the 
projected distributions F(xi), G{x\) or F(x 2 ), G(x 2 ) are not (F(xi t2 )\G(xi t2 )) ^ 0. 
When we project f{x\,x 2 ), g(x\,x 2 ) on the effective orthogonalization variable x, 
defined by the direction between the distribution "centers" , the resulting projected 
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V-''' g(Xi,X 2 ) 



F(xi) G(xi) Xl ^> 

Figure 1: Scheme of the distributions f(xi, x%), g{x\, x%) for two species, quasi-orthogonal ((f\g) ~ 
0) in the parameters x\, x 2 - The projected distributions F(x\), G{x\) or F(x 2 ), G(x 2 ) overlap 
{F{ x H2))\G{xn2))) 7^ 0, while the distributions F(x), G{x) resulting by the projection along the 
effective orthogonalization variable x are still nearly orthogonal (F(x)\G(x)} ~ 0. 

distributions F(x), G(x) become also nearly orthogonal (F(x)\G(x)) ~ 0. This 
reveals the OP-built correlations between parameters x±, x 2 in the two species. 



Orthogonalization in sympatric speciation 

Investigation of OP in species is a demanding task, due to the formation of the 
complex correlations between species parameters and the different involved levels. 
OMAL processes can be often observed in evolutionary or functionally close species. 
Very interesting from the point of OP is the sympatric speciation, with its many 
mutually interconnected aspects, like the presence of competition or the role of 
assortativity in sexual interbreeding (Doebeli (1996b), Kondrashov & Kondrashov 
(1999), Dieckman & Doebeli (1999)). The activity of OMAL could help us to under- 
stand some problems of sympatry, like its possible existence on a "flat landscape" 
solely due to competition (Rosenzweig (1978)). 

SPECIATION EQUATION 

Here, we model sympatric speciation on the macro-level, where we investigate the 
presence of OMAL. For simplicity, we do not explicitly model here the parallel 
fixation on the micro-level. Only sexually reproducing organisms are considered, 
but the results could be extended also to clonal species. We describe the total 
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character distribution /(x) by the new equation 




= ft($,x) J J dsn dx u Sfaxn) (/(x-xj) /(x + x /7 )) a -^g| , (4) 



and assume that /(x) splits into several sub-populations /i(x) for individual species 
with a limited interbreeding. The front term on the r.h.s., representing growth of 
the species, is controlled by the effective resources 1Z, with abiotic and biotic compo- 
nents, and dependence on factors from the set $, like temperature or radiation. It is 
also influenced by the sexual function S with a second order mating (two partners) 
in a power a. The last term in Eqn.(f|) represents dying of individuals with a time 
t(x), due to natural reasons, predator-prey coupling or catastrophes. It can also 
reflect possible extinctions, since /(x) would keep a constant norm for r(x) ~ oo. 



The resources function 1Z plays the role of a "fitness landscape" , considered in other 
speciation studies (Peliti (2000)). We use it in the form 



where competition can be largely varied. Here, TZq(<&, x) represents all the resources 
available for an individual with parameter x, in the absence of other individuals, 
and the external parameters We simply assume that it is given by the Gaussian 
function 7?-o( x ) = Ro e~' x_Xfl ' 2 ^°" fl - )2 . Exploitation of the resources is controlled by 
the exponential function in Eqn.([|) with the distribution /(x), which is competition- 
broadened by the function 1Z C . This broadening expresses the fact that individuals 
with parameter x' can eat the food already consumed by individuals with parameter 
x. If we assume that TZ C is also Gaussian, its width oq determines the effective "dis- 
tance" of the consumption/competition, and additional nonlinearity of this process 
can be tuned by the parameter f3. 

Roughening of the landscape lZo(&, x) reflects the discrete character of species, 
which form the resources, and the structured ability of other species to consume 
them. In principle, both can be related to a "limited scaling" in the species pa- 
rameters. For example, individuals of some species cannot get just smaller without 
qualitatively changing their properties, like the ability to hide by digging a hole, be- 
cause of the increasing danger from predators, as they become smaller. In this way, 
complex correlations in species parameters are built by OMAL, which are roughen- 
ing the landscape TZo(<&,x). We do not explicitly consider here the food-web struc- 
ture, where landscape roughening could propagate, and neglect also predator-prey 
dynamics. 



RESOURCES FUNCTION 



ft($,x)=72o($,x) exp{-(/ rfx ft c ($, x - x)/(x)) /3 } , 



(5) 
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SEXUAL TERM 



In the growth term from Eqn.(^j), we assume that the progeny population /(x) 
receives the properties x from the parents with distributions /(x — xj), /(x + x/j). 
In the mean heritage approximation, only parents with x/ = kjj contribute to this 
solution. The dependence of the birth rate on the distance of parents parameters 
X = x/ + x// is given by the sexual correlation function S (X, £) . Declination 
of the progeny properties from this mean solution is described by the difference 
£ = x/ — x/7, which originates in mutations and other effects, like gene mixing and 
activation. We can again assume that S has a Gaussian form 

Ixiv4-KI 2 A^ 

where the width of sexual interests and the declination om can be partly tuned 
by OMAL. The factor a in Eqn.(|j) reflects the "power of interbreeding". Speciation 
is rather sensitive to the departure of the equation from the quasi-linearity a ~ 0.5, 
as we discuss below. The nonlinearity could be also incorporated in a more general 
version of the decay term — /(x)/r(x) in Eqn.(Q), where it would reflect a decreased 
defense ability of small size populations. 

Numerical Studies 

In numerical studies, we use projected distributions with just one scalar property x 
for the effective orthogonalization variable (see Fig. [I|), in which the species overlap 
the least. We assume, for simplicity, that the projected distribution F(x) (called 
here also f(x)) follows the same Eqn. f| as the full distribution f(xi,X2,...,x n ). 
We define the formed species i as separate sub-populations fi(x) in the character 
distribution f(x), similarly as in (Doebeli (1996b), Drossel & Mac Kane (2000)). 
The species i, j are considered to be orthogonal if their final steady-state distribu- 
tions fi(x), fj(x) do not overlap (fi\fj) ~ 0. The sub-populations of the consump- 
tion/competition function fc(x) can still overlap (fa\fcj) 0; which leads to the 
fact that the species separation in the variable x is mostly controlled by the finite 
consumption/competition width ac- 

We transform in Eqn. ^ the internal variables (xi, xn) or equivalently (X, £) 
to the new variables (z, z), where X = 2z, £ = 2(x — z). Then the equation reads 

*r - / / » * - "^P 1 & - * * - ^ . m 

where, the resources function is 



7Z(x) = R e~ x t«R-°oVoW) ; f c ( x ) = / dx = f(x) 
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In Eqns. (0-|8|) we use the values S° = 0.2 and ac = 4 to easily present numerical 
results for different functions. 

PRACTICAL EXAMPLES 

We use Eqns.(0-||) to study under which conditions sympatric speciation can be 
obtained on broad Gaussian resources, and pay also attention to the role of their 
local inhomogeneities. We also specify the role of OMAL in the obtained solutions. 

In Fig.|2| we present the steady-state (globally stable) numerical solutions of 
Eqns.(0-||) for species formed by sympatric speciation. The dashed and dash-dotted 
curves correspond to the free and unused resources without and with the competition 
exponential from (§), respectively. The solid and dotted curves represent the total 
population f(x) and the consumption/competition function fc(x), respectively. The 
situations with one to four species are obtained for the Gaussian resources of the 
width or = 0.2 and strengths Ro = 0.5, 1, 2, 4. The consumption/competition width 
is ac = 0.1, the sexual and mutation widths are chosen the same as = au = 0.015, 
the power parameters are a = j3 = 0.6 and the lifetime is r = 100. 

The broad Gaussian resources can be seen as "quasi-flat", since a^ is larger 
than other parameters in the model. When the richness of resources Ro is rather 
small, only ac is close in size to or. As Ro increases, an effective width of resources 
Tir (on which species can survive) becomes larger than or, so the resources can 
feed more and more species. If an alone (not Ro) is enlarged to get a nearly flat 
landscape, the number of species grows, but their separation is always given by ac- 
If ac does not vary much, the resulting self-organized speciation quasi-periodicity 
as a function of x determines a universal speciation volume. In the present broad 
Gaussian landscape, with an exponential form of competition (see (|)), speciation 
can be obtained only if the power of interbreeding is slightly nonlinear a > 0.5. This 
adds to the speciation conditions on a quasi-flat landscape. 

The unused resources (— . — .) in Fig|| are flattened by the total distribution 

fc{x), which largely copies the shape of the resources ( — ), even though 

the total distribution f(x) has sharply separated sub-populations. This shows an 
efficient sharing of resources by the orthogonalized consumption needs of the species 
(Roughgarden (1976)). We have used Eqn. |3]to calculate the relative overlaps r C ij = 
(fa\fcj)/(fc\fc) of the individual contributions fa{x) to f c (x). These fa(x) can 
be obtained from Eqn. ||, where we substitute the individual distributions f(x). 
The situations with n = 2 — 4 species in Fig.|2| give: n = 2, r\2 = 0.296; n = 3, 
ri 2 = 0.219, = 0.045 and n = 4, r 12 = 0.17, r 13 = 0.062, r 14 = 0.008, r 23 = 0.213. 
The values rcij show the steady-state strength of coupling of the orthogonal species 
{(fi\fj) ~ 0). The related pressures can be relaxed by inhomogeneities in the broad 
Gaussian resources. 
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Figure 2: Species population / (competitive population fc) in a sympatric speciation is shown by 
full (dotted) lines as a function of the effective parameter x. The dashed (dash-dotted) lines rep- 
resent the free (competitive) landscape of resources. The plots with one to four species correspond 
to strength of resources Rq = 0.5, 1, 2, 4 and other parameters in the text. 
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Note also that the final distributions fi(x) in Fig.|] do not depend on the size 
of the initial Gaussian population used in iterations of Eqns.(0-^), unless this is very- 
small. Then, for the present parameters, initial populations which are by 2-3 orders 
smaller than the final population tend to follow the single-peaked landscape. This 
means that the final solution gives 2n — 1 species where 2n species would normally 
appear, so the initial condition is "frozen in" the final state of a reasonable stability. 
For even smaller initial populations species die out, due to finite lifetime r. 

Stepwise speciation resembling our results has been observed, for example, 
in stickleback species in small lakes of the coastal British Columbia (Schluter & 
Mac Phail (1992)). Speciation driven by the increased strength of resources could 
be also the reason of diversity gradients of marine life, as induced by the solar 
radiation in coastal regions (Roy et al. (1998)). Similarly, richness gradients of other 
species observed towards the equator (Rapoport (1975)) could be the consequence 
of enlarged radiation, temperature and food zones. 
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Figure 3: (top plot) The maxima % of the populations, for the species with parameters in Fig. ^, 
as a function of the resources strength Rq. (bottom plot) The same dependence for the maxima 
of the species populations Jm, with the number of peaks denoted. 

In Fig.|3| we show the size fu and position i M of the sub-population maxima 
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as a function of the resources strength Rq and other parameters as in Fig. |^. In 
the top plot for xm, we can see that after a new species emerges the positions xm 
for the peaks slightly move. Since the Gaussian resources rapidly fall down at the 
tails, the effective width increases very slowly with Rq, and so does the number 
of species n for n > 3. For example, the fifth species appears only at Rq m 13. In 
the bottom plot, we show the dependence of the maxima /m on Rq, and add the 
number of peaks in each curve. As Rq grows and the number of species increases 
from odd to even, the main maxima jump in values for this shape of resources, while 
in transitions from even to odd the changes are reasonably continuous. The separa- 
tion of species parameters (character shift) shows that OMAL is accompanied by a 
pressure between the coexisting species (Slatkin (1980)), induced by the overlap of 
their fa(%) distributions. The activity is especially large close to speciation events, 
where the species parameters move. Experimentally observed pressure induced vari- 
ation of the character shift, as the species move between an island and a continent, 
is discussed for example by Schluter (1988). 

Population branching could be also studied in a time region (Kisdi (1999)). In 
our work, globally stable solutions in Fig.Q are calculated by using at each point R 
a Gaussian initial population comparable in size with the final value. We will also 
study a time evolution, but limit to the investigation of speciation for an adiabatic 
time evolution of the resources richness Ro; we slightly change Ro in each time run, 
and start from the populations of the previous Ro. The results for the species posi- 
tions shown in Fig.^, where thin doted lines correspond to the solutions from 
FigfJ. In the upper plot, we show by thick solid (dashed) lines runs with increasing 
R starting from the situation in the single-species (triple-species) region. In the 
former case, the first speciation occurs delayed in the double-species region, while 
three species would appear only deeply in the region where four species should be; a 
start from the double-species region follows this solution as well. If we start from the 
triple-species region, the number of species keeps unchanged for the used parameters. 
If the resources become poorer, the situation is rather different, as we show on the 
bottom plot. A start from the four-species region gives evolution skipping over so- 
lutions with three species, but it "smoothly transfers in advance" to a solution with 
two species, and the solution with one species appears delayed. Starting from the 
triple-species region, the evolution skips over solutions with two species, and later 
jumps to a solution with one species. Thus, if Ro increases, all numbers of species 
appear, but they are delayed, while if Ro decreases, sequences with odd or even 
number of species seem to be realized. Since the frozen nonequilibrium solutions are 
locally stable with respect to perturbations, the system experiences hysteresis in the 
number of species if R slowly oscillates with a large enough amplitude. It is likely 
that the species parameters would change by OMAL, to partially follow the slow 
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evolution of resources and bypass the mounting pressures from a fixed number of 
species. The crucial aspect is the ratio of timescales of these effects. In nature, such 
a behavior can be also largely stabilized by year seasons and other random effects, 
but if we take into account the evolution on micro-level one could essentially face 
only irreversible behavior. 
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Figure 4: Positions xm of the species as the resources richness Ro is slowly varied. The upper 
(lower) plot shows situations for increasing (decreasing) R . Thin doted lines correspond to the 
globally stable solutions from Fig.|| In the upper plot thick solid (dashed) lines present runs 
starting from regions with a single (triple) species. In the lower plot the solid (dashed) lines 
present runs starting from regions with a triple (quadruple) species. In all cases speciation freezing, 
dependent on history, can be observed. 

In Fig. | we show speciation in the previous landscapes, modified by additional 
disturbances. The disturbances can trigger speciation at smaller competition pres- 
sures, induced by the fa distributions; speciation could even take place for a = 0.5, 
which would not be possible on the broad Gaussian landscape. In the upper plot, 
we take the situation with one species in Fig.§, and subtract from the resources 

2/2 2/2 

Ro e~ x > a R the exponential R\ e~ x ' a p with R\ = 0.1, o~p = 0.07. The resulting 
landscape flattening can split the single species in two sub-populations, since the 
resources become more efficiently explored by two species. We have also found such 
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Figure 5: The situations from Fig.|| presented for axperturbed Gaussian landscape. The top plot 
shows how the single species from Figj2] breaks in two, when the resources are flattened. Next, 
one of the side species, from the third plot in Fig.§, shifts and decreases in size to fit the new 
landscape. In the third plot, for the same parameters, but the sexuality width as — 0.035, the 
side peaks dissapear, but the resources are not well explored. This is improved on the last plot, 
where, instead of as, the consumption width is broadened to the value ac = 0.15. 



14 



splitting/joining jumps for the case of two or three species. For R = 1.5 the num- 
ber of species goes from three to two, while for R = 2.5 it goes from three to four, 
copying thereby in both cases the flatter shape of resources. 

From a theoretical point of view, we could also consider that x represents one of 
the parameters $, like temperature, radiation or humidity. Then, the resources show 
how much food for the species is available at this parameter x, and other changes are 
analogous. Therefore, we should in principle be able to describe certain situations 
with allopatric speciation. For example, splitting of a Drosophila population, similar 
to the top part of Fig.|5|, was observed on a microsite in a valley with two differently 
irradiated slopes, giving a radiation/temperature space gradient corresponding to 
our parameter x (Harry et al. (1999)). A sharply different character of the slope 
sides could give flatter resources in the middle radiation/temperatures with side 
maxima as in Fig.[5[ 

In the second plot of Fig.|5|, we present the population for the situation with 
three species in Fig.|2|, where we subtract, similarly as before, the shifted exponential 
Ri e~( x+Xl ^ 2 / a F with Ri — 1, X\ — 0.1 and o> = 0.07. Due to this asymmetry of the 
resources, the first side species largely shifts and decreases, in order to efficiently 
explore the new landscape, but the total number of species still remains the same. 
This inhomogeneity-shifted species competes a little with the other species. In the 
next plot, we use the same resources, but we broaden the sexual width to as = 0.035. 
This decreased assortativity fully prevents speciation, and the resources are poorly 
used in the tails. Much better exploration is obtained in the last plot, where, 
instead, the consumption width is increased by 50 % to oc — 0.15, so the single 
species population is also larger. The width om is also crucial, since it can broaden 
and smooth the population. For example, the single peaked distribution in the third 
plot of Fig. |^ can be also obtained if om is increased to the value 0.035. 

GENERAL TUNING OF PARAMETERS 

Other parameters in Eqns. (0-|D could also interplay in a sensitive way. For example, 
the power a can fluctuate around the value a ~ 0.5 to promote/suppress speciation. 
The power (3 can reflect the way of food consumption as well as possible deterrence 
actions aimed in protecting it (screening). For (3 < 0.5 (/3 > 0.5) side populations 
become suppressed (enhanced) with respect to the middle populations. The absolute 
strength of consumption/competition ac can be tuned to control the size of the 
population. Finally, the time r maintains the speciation pressure and keeps the 
population stability (Smith & Brown (1986)). The term with r in Eqn.(f|) diminishes 
individuals present in regions between the landscape maxima. It also gives sharper 
and a more regular speciation. Without competition and with large r, the population 
would simply follow the landscape, while with competition, it would not speciate 
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sharply until r is decreased. 

Recent studies point to the increased role of food-web resources in species 
survival (Amaral & Meyer (1999)). Proportionality between the number of species 
and the resources richness, obtained above, can be also supported in different ways. 
For example, in poor resources, generalist species might be triggered by broadening 
ac and as- Clusters of species can be extinct, and surviving species can broaden 
their distributions, so that correlations formed by OM AL become decreased. Since 
mutations are less suppressed, speciation could eventually be also promoted in these 
conditions. If the resources get richer, individualists could do better. Species can 
narrow their consumption needs ac-, but they need to live efficiently in other ways. 
Consequently, they also narrow their sexual interests as and trigger further specia- 
tion. Since the populations largely copy resources, their richness does not make life 
easier for the numerous individuals. This drives the lasting self-organized evolution. 
OMAL can tune all these parameters in a limited amount to balance coexistence 
of species with an efficiency exploration of resources. The timescale of this activity 
goes over many generations, when these changes get fixed on the micro-level. 

FIXATION OF OMAL 

We noted that neutral sections of proteins, with a relatively stable composition 
in each individual species, might be used for transcription of OMAL. This can 
support selectionist views on the problem of molecular evolution clock (Ohta (1996)), 
which underlines the quasi-deterministic origin of observed constant mutation rates 
in neutral sections of proteins. Recent observations show, that when bacteria are 
placed in a new environment, some of their genes start to mutate faster (Sniegowski 
et al. (1997), Schmid Sz Tautz (1997)). Growth of RNA in vitro also shows that 
under strong selections the molecules go to the same limit form (Strunk & Ederhof 
(1997)). These experiments point to a possibly deterministic role of the environment 
on the proteins, where OMAL takes an active role. 

We can imagine, that tuning of protein catalytic strength can be regulated by 
other proteins, enforced by the complex molecular environment, evolving under the 
pressure of globally changing conditions on the macro-level. If these possibilities are 
exploited, and a larger change is needed, a jump in neutral positions of the tuned 
proteins could be induced by OMAL, which is analogous to avalanche ticking in 
other systems (Wu et al. (1993)). This ticking could largely go in parallel with the 
formation of new species, shown for example in Fig. j|. 

These ideas could be tested on bacteria, with short reproduction periods. A 
working scenario could be obtained by modifying some neutral protein region, and 
keeping the remaining gene-pool frozen by a periodical exchange. If OMAL is rele- 
vant, the protein would eventually evolve in many bacteria generations to the same 
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or a similar form. The problem is if OMAL can be reached in artificial laboratory 
conditions. Such experiments could also reveal to which extent one can speak about 
proteins in nonequilibrium with the cell environment. 

Summary 

We have introduced the concept of species orthogonalization on the micro and macro- 
level. OMAL is related with the need for formation of complementary functions 
and habits of different biological species living in close contacts. OMIL separates 
functions of molecular species, and assists OMAL on the micro-level. We have 
shortly discussed the ways of evaluating orthogonalization processes. 

As a practical example, we have investigated OMAL in sympatric speciation, 
but limited so far our study to the macro-level. To this goal, we have found a new 
population equation, and presented it in several numerical examples. Its solutions, 
for the projected distribution with one effective variable x, show that the number 
of species grows in richer resources, divided between the species. This could help 
to explain the increase of species diversity close to the equator, where sun radiation 
provides richer resources. We have also observed speciation dependent on history, 
which might be useful in modeling of ecology, and discussed possible ways of OMAL 
fixation on the micro-level. We expect that orthogonalization mechanisms are active 
in other systems with competitive structures. 
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